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ABSTRACT 

We report ^100 pc (0'.'3) resolution observations of (sub) millimeter HCO + and CO lines in the 
ultraluminous infrared galaxy Arp 220. The lines peak at two merger nuclei, with HCO + being more 
spatially concentrated than CO. Asymmetric line profiles with blueshifted absorption and redshiftcd 
emission are discovered in HCO + (3-2) and (4-3) toward the two nuclei and in CO (3-2) toward one 
nucleus. We suggest that these P-Cygni profiles are due to ~100 km s _1 outward motion of molecular 
gas from the nuclei. This gas is most likely outflowing from the inner regions of the two nuclear disks 
rotating around individual nuclei, clearing the shroud around the luminosity sources there. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: individual (Arp 220) 
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1. INTRODUCTION 
The two merg er nuclei of Arp 220, about 300 pc apart 



r p 220 t 

on the sky (see IScoville et ail fl998. 1"=361 pc), gen- 
erate at least 40% of the total luminosity of the ul- 
traluminous infrared galaxy (Ls-iooo^m = lO 12 ' 2 ^©; 
iSoifer et al.lll999l ). Each has a rotating gas disk of about 
100 pc extent, the rotation axes of which are misaligned 
with one another (Sakamot o et alj fl999). The central 
50-80 pc of the brighter western nucleus has a bolometric 
luminosity of > 2 x 1O 11 L accord ing to sub-arcsec imag- 
ing of submillimeter continuum (|Sakamoto et all [20081 
hereafter IS08I ). There is als o a kpc-scale molec ular disk 
surrounding the two nu clei flScoville et alJl!997f). IR ab- 
sorption analyse s (e.g.. iDudlev &z Wvnn- Williams! [1997T . 
lHaas et aLll2001h suggest the dominant luminosity source 
of Arp 220 to be compact and deeply buried, as ex- 
pected in the two nuclear disks, rather than widely 
spread, as expected in the outer disk. Both nuclear disks 
host starbursts wi th dozens of radio supernova remnants 



Parra et a l. 2007) and may also contain a buried quasar 



Soif er et al1ll98i [D ownes fc Eckart 2007 and references 
therein; |S08|). 

Properties of the nuclear disks, including their struc- 
ture, kinematics, and gas conditions, must be elucidated 
in order to understand fully the nature and evolution of 
activity in Arp 220. In this Letter we report our high- 
resolution observations of Arp 220 in (sub) millimeter 
molecular lines, including its first subarcsecond resolu- 
tion imaging in HCO + and the first detection of P-Cygni 
line profiles toward the two nuclei. 
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2. SUBMILLIMETER ARRAY OBSERVATIONS 

We observed CO(3-2) and HCO+(4-3) simultaneously 
in 2007 May and HCO+(3-2) in 2008 August using the 
most e xtended co nfiguration of the Submillimeter Array 
(SMA; IHo et al.ll200l 8 . We also used the most com- 
pact configuration in 2009 January to measure the to- 
tal HCO + (3-2) flux. System gain was calibrated using 
two nearby quasars, J1635+381 and J1613+342. The 
passband over our 2 GHz-wide sidebands was calibrated 
using other brighter quas ars. The data were reduced 
mostly the same way as in lS08l Continuum in each side- 
band was collected from the channels outside of radial 
velocities between 4800 and 5900 km s" 1 in CO(3-2), 
HCO+(4-3) and (3-2), and HCN(4-3) lines. The con- 
tinuum was then subtracted in visibilities from line data. 
Phase self-calibration was made using the continuum of 
Arp 220. Line data were binned to 30 km s _1 resolution 
and have velocities in the radio convention with respect 
to the LSR. Our spatial resolution is about O^ or 100 
pc and our CO(3-2) data are twice as sensitive as those 
in IS081 Since CO(3-2) was placed at the band center, 
our velocity coverage for HCO + (4-3) in the image band 
is limited to 5180 - 6840 km s -1 . 

3. RESULTS 

3.1. Compact HCO + and Extended CO 

We detected HCO + , CO, and A~l mm continuum 
emission associated with the two nuclei (Fig. []}. The 
western nucleus is more luminous than the eastern, by a 
factor of 2-3 in HCO + and continuum. The concentra- 
tion of HCO + emission to the two nuclei is real because 
we detected most of the total line flux in Arp 220 (Table 
Q}. This is in contrast to CO (3-2) where only 23% of the 
single-dish flux was recovered even though CO (3-2) and 
HCO + (4-3) had almost the same u-v coverage. CO (3-2) 
emis sion is indeed extended on scales greater than 2" in 
[S08l . Thus HCO+(4-3) and (3-2) emission of Arp 220 is 
mostly from the two nuclear disks while the majority of 

8 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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CO (3-2) comes from the outer disk around them. More- 
over, in each nuclear disk, the HCO + emitting region is 
more compact than the CO emitting region judging from 
their half-peak sizes in Fig. [T] 

The peak brightness temperatures (Tb) of the lines at 
or around each nucleus, listed in Table [TJ indicate warm 
molecular gas in the nuclear disks. The CO peaks are 
about 50 K and HCO+ about 20 K. These high T b and 
high HCO + -to-CO intensity ratios in the nuclear disks 
imply that molecular gas exists there in conditions fa- 
vorable for the transitions of high critical density. 

3.2. P-Cygni Line Profiles at Both Nuclei 

The line spectra at the two nuclei, in Fig. [21 show both 
emission and absorption. In particular, all the observed 
lines except CO (3-2) toward the western nucleus show 
blucshiftcd absorption and redshifted emission, i.e., P- 
Cygni profiles. This feature is absent in our calibrator 
spectra. 

The western nucleus has a narrow, deep absorption 
feature at 5340 km s -1 and emission at higher and lower 
velocities in our CO data. This profile is consistent with 
the previous CO(2-l) profile in lDownes fc Eckartl ([20071 ) 
except that the absorption is deeper in our CO(3-2) data. 
The deep absorption is almost at the systemic velocity 
of the nucleus, marked with a dashed line in Fig. [21 
as expected if the CO-absorbing gas is in circular rota- 
tion in the western nuclear disk. The systemic velocity, 
T4y S (W) = 5355±15 km s _1 , was estimated from the cen- 
tral vel ocity of the CO profile in the low-resolution (0'.'5) 
data of IS08L measured at half- maximum. The HCO + (4- 
3) and (3-2) lines are also in absorption at the systemic 
velocity. In addition, both HCO + lines have absorption 
at blueshifted velocities while most of their emission ap- 
pears at redshifted velocities. The HCO + (3-2) profile 
shows multiple velocity components in absorption; the 
minima at 4890, 5220, and 5340 km s _1 are at 3.6er 4.6er, 
and 4.6(7, respectively. Absorption deeper than 3<r is seen 
down to Ky S (W) - 495 km s" 1 . 

The eastern nucleus has very asymmetric line profiles 
with respect to its systemic velocity, which we estimate to 
be T4 ys (E) = 5415±15 km s _1 in the same way as for the 
western nucleus. In all the observed lines, absorption is 
predominantly seen at or blueward of the systemic veloc- 
ity while most emission is in redshifted velocities, making 
the profiles P Cyg type. The deepest absorption is at 45- 
75 km s _1 below V sys (E) and is 2.8ct, 4.5ct, and 3. Oct in 
HCO+(4-3), (3-2), and CO (3-2), respectively. 

3.3. Absorption Depth 

Absorption (i.e., negative intensity) in our continuum- 
subtracted data means absorption of the continuum, be- 
cause line self-absorption alone cannot produce negative 
intensity. In order to absorb continuum the excitation 
temperature of the foreground gas must be below the 
brightness temperature of the background continuum. 
The (deconvolved) western nucleus has Tb= 9 0-16 K 
and a size of 50-80 pc in 860 /zm continuum (|S08l; fo r 
the SED of the two nuclei see [Matsushi ta et al.l 120091) . 
The high Tb, due to high dust opacity in submillimc- 
ter, makes line absorption against dust continuum more 
likely to occur (and deeper) than at longer wavelengths. 

The maximum apparent optical depth of each absorp- 
tion line is in the range of 0.4-1.0 except the HCO + (3-2) 



optical depth of 3.7 (or >1.4 when allowing for 1 a error) 
toward the eastern nucleus. The total line absorption 
/ r a dV integrated over the velocities where r a > 0.1 is, in 
the order of HCO+(4-3), HCO+(3-2), and CO(3-2), 75, 
150, and 50 km s toward the eastern nucleus and 85, 
82, and 15 km s _1 toward the western. The absorption 
profiles vary among the lines, more notably toward the 
western nucleus. There are two caveats to keep in mind. 
First, our observations provide lower limits to total ab- 
sorption because emission from the absorbing gas itself 
and from ambient (i.e., non-absorbing) gas in our beam 
masks the absorption. Second, our synthesized beam is 
larger for HCO+(3-2) than for HCO+(4-3) and CO(3- 

The absorbing column density N^ 2 is on the order of 
10 22 cm -2 for the CO absorption and 10 23 cm -2 for the 
HCO + if the absorbing gas is in 50 K LTE and abun- 
dances are [CO/H a ] « 10" 4 and [HCO+/H 2 ] « lO" 8 . 
The high abundance adopted for HCO + is consistent 
with recent PDR and XDR models (Meiicrin k fc Spaansl 
120051) . as well as models of HCO + enhanceme nt in pro- 
tostellar outflow shocks (|Rawlings et al.l l2004). The ab- 
sorbing gas will have a depth on the order of 1-10 pc if 
its density is at the CO (3-2) critical density of 10 4 cm~ 3 . 
The actual absorbing material can be distributed over a 
longer length if it is inhomogeneous as in numerical hy- 
drodynamical sim ulations of multi-phase gas obscuring 
active nuclei (e.g. JWada &: Norman! [2002f) . The column 
densities of the absorbing gas calculated here have large 
uncertainties because the absorption data provide lower 
limits, gas excitation conditions are poorly known, and 
molecular abundances can significantly vary. Still, it is 
possible for the absorbing column density to take these 
values because they are smaller t han the total column 
density of the nuclei (~10 25 cm~ 2 ; IS08T ). The remarkable 
detection of vibrationally e xcited HCN in absorption in 
Arp 220 (|Salter et al.ll2008l ) suggests that the excitation 
of HCO + may also be more unusual than described by 
our simple analysis. More realistic models of the excita- 
tion and radiative transfer of both molecules, including 
radiative coupling to the intense infrared continuum, will 
be explored in a future publication. 

3.4. Gas Motion — Rotation and Outflow 

Rotation of gas around each nucleus is evident in 
the CO(3-2) position-velocity diagrams (Fig. [3]) as the 
negative to positive shift of V ~ V sys across the nu- 
cleus. The CO emission ranges from 5000 to 5900 
km s _1 around the western nucleus. The high-velocity 
gas around 5800-5900 km s _1 is detected at >3 a for 
the first time. The fall-off of the rotation velocity away 
fro m the nucleus agr e es wit h the CO (2-1) observations 
by iDownes h Eckartl ([2007) and suggests truncation of 
the mass distribution of the nucleus beyond a <100 pc 
core. 

Radial motion of gas, in addition to the rotation, is sug- 
gested by the blueshifted absorption mentioned above. 
The absorbing gas must be in front of the continuum 
core of each nucleus and moving away from it toward us. 
Although absorption only tells us the gas motion along 
our lines of sight to the continuum nuclei, the detection of 
similar blueshifted absorption on both nuclei whose nu- 
clear disks are misaligned suggests that the gas motion is 
more likely a radial outflow in most directions than non- 
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circular motion only along our sight lines. The redshiftcd 
emission in the P-Cygni profiles must then contain emis- 
sion from gas on the far side of the continuum nuclei and 
moving away from the nuclei and us. The other emis- 
sion in the spectra must be from the gas rotating in the 
nuclear disks. 

The typical outflow velocity of the absorbing gas is 
probably about 100 km s -1 because the deepest absorp- 
tion in the P-Cygni spectra is between (V sys — 45) and 
(Vgyg — 135) km s _1 . High-velocity components up to 
^500 km s _1 are suggested toward the western nucleus 
by the HCO + (3-2) absorption with large velocity offsets. 
The absorbing gas is turbulent or has a velocity gradi- 
ent across each continuum nucleus because absorption is 
also seen at slightly positive velocities with respect to 
systemic. The mass of the outflowing gas and the mass 
outflow rate would be on the order of 5 x 10 7 M Q and 
100 Mq yr _1 if the outflowing gas with a column den- 
sity of 10 23 H2 cm~ 2 is mostly within a radius of 50 pc 
(see §4.1|) and has an isotropic outward velocity of 100 
km s _1 . Note however that these numbers strongly de- 
pend on the geometry of the gas flow. 

4. DISCUSSION 
4.1. Location of the Absorber 

The blueshifted absorbing gas is probably located in 
the two nuclear disks. This is likely because HCO + emis- 
sion is concentrated to the two nuclear disks. The nuclear 
disks thus have significant column densities of HCO + ex- 
cited to the J=3 and 4 levels, whether the excitation 
is collisional (with critical densities of 10 6 -10 7 cm~ 3 at 
50 K) or by radiation. The J=3 HCO+ in front of a 
bright source causes HCO + (4<— 3) absorption. Since col- 
umn densities of excited HCO + must be less outside the 
nuclear disks according to our observations, any external 
gas is unlikely to be the dominant absorber. 

In our preferred model, each nuclear gas disk has an in- 
ner region with significant outflow motion and an outer 
region dominated by rotation. This configuration can 
partly explain why the P-Cygni profiles are more promi- 
nent in HCO + than in CO. Excited HCO + is more con- 
centrated toward the outflowing disk centers, while for 
extended CO more emission from the rotating (i.e., non- 
outflowing) disk likely masks absorption by the outflow- 
ing gas. It is also likely that the lower excitation tem- 
perature of HCO + expected from its larger Einstein A 
coefficients makes it a better absorber than the CO at 
the same location with a higher excitation temperature. 
The rotation-dominated part of the nuclear disk is likely 
the si te of possible weak HNC maser that lAalto et al.l 
(2009) suggested for the western nucleus on the basis of 
a narrow and bright spectral feature at the systemic ve- 
locity. 

A less likely location of the absorbing gas is a few 100 
pc or more from the two nuclear disks. The gas may 
be a part of the merger's large-scale (> 10 kpc) super- 
wind perpendicular to the kpc-scale outer mo lecular disk 
(jHeckman et al.lll987l : iMcDowell et alJl2003l ). The 5340 
km s _1 absorption toward both nuclei could be explained 
by this common envelope model. However, the total mass 
of such a gas structure would be quite large; ~10 9 M Q 
for an expanding envelope of 500 pc radius and a col- 
umn density of _/Vh 2 ~ 10 22 5 cm~ 2 . Such an expanding 



molecular envelope has not been detected in emission. 

4.2. Relation to Previous Observations 

Observational signs of outflows from individual nuclei 
in Arp 220 have been scarce despite ample observations 
of the merger's superwind at larger scales. A bicon- 
ical distribution of OH maser around the western nu- 
cleus alo ng its rotation axis (i .e., N-S directi on) was ob- 
served bv lRovilos et all (|2003l their Fig. 4). lAalto et all 
(2009) suggested that this emission as well as the faint 
HNC emission they found to be extended in the same 
direction may be due to an outflow from the nucleus. 
While kinematical evidence for the outflow is absent 
or incon clusive in these su barcsecond resolution obser- 
vations, iBaan et al.l (|1989h detected highly blueshifted 
OH maser emission with a ~2' beam and attributed it 
to an outflow at 10-100 pc from the luminosity source 
with wind acceleratio n modeling. In H I observations by 
iMundell et"ai1 (|2001h . velocities of the main absorption 
toward the nuclei are V(E, main H I) = 5396 km s _1 and 
V(W, main H I) = 5325 km s _1 after conversion from 
optical-heliocentric to radio-LSR. They are marginally 
blueshifted, by 19 and 30 km s _1 , from our CO-based 
systemic velocities. Also, the H I absorption toward 
the western nucleus has a subcomponent at about —80 
km s _1 from the main component. These may be related 
to the molecular outflows that we suggest. 

4.3. Nuclear Winds and the Nuclear Disks 

The P-Cygni profiles toward the two nuclei provide a 
new clue to the energetic nucleus of Arp 220. The nuclear 
disks need to be reinvestigated in light of the suggested 
nuclear winds. 

The compact and luminous sources in the nuclei prob- 
ably make the outflows inevitable. For example, super- 
novae inject 3 x 10 50 J of mechanical energy into the 
western n ucleus in 1 Myr for th e supernova rate of ~3 
yr _1 that lLonsdale et all (|2006h observed. (The rate is 
^1 yr _1 for the eastern nucleus.) If 10% of this is passed 
to kinetic energy of 10 8 Mq of gas, the gas velocity will 
be 500 km s _1 . Gas at shallow parts of the gravitational 
potential or getting more energy than average will es- 
cape the nucleus. Although the outflow may be mostly 
ionized, it will entrain molecular gas on its way. Radia- 
tion pressure on dust also drives an outflow in the case 
of a very compact starburst or an AGN buried in dus t 
(|Scoville et alil2001l: IScovfli7ll200l iMurrav et al.ll2005T) . 
The momentum of a 5 x 10 7 Mq gas outflowing at 100 
km s _1 is equivalent to the total momentum of photons 
that a 0.5 x 10 12 L Q source radiates in 0.5 Myr. Such an 
outflow may be more isotropic near its origin than the 
outflow caused by a disk starburst and would be more 
significant in the inner region of each nuclear disk be- 
cause the radiation pressure declines with radius as r -2 . 

The winds are important for the evolution of the nu- 
clear disks because they remove the shroud of and fuel 
for the buried energy sources. The small crossing time of 
the nuclear winds over each disk, <1 Myr, is suggestive 
in terms of the evolution. It is either that we are seeing 
very young winds or that the gas motion in each disk is 
not purely outward in every direction despite the large 
opening angles of the outflows inferred in £ 13.41 For the 
latter, it is conceivable that high-density gas in the mid- 
plane of each nuclear disk is little affected by the outflow 
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and even has inward mot ion due to viscosit y to re plenish 
gas to the center, as in iWada fc Nor"rnanl ([2001 . The 
evolution of the nuclear gas disks is then determined by 
the balance of the gas accretion, gas outflow, and gas 
consumption by star formation and AGN until the coa- 
lescence of the nuclei and of the gas disks. 
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Fig. 1. — Integrated intensity maps of HCO + (4-3), (3—2), and CO(3— 2) lines and 0.87 mm continuum. The line moment-0 maps were 
made using positive signals only. The lowest contour is at 2cr and the n-th contour is at 9.5n, 4.3n, and 10. 3n Jy beam - 1 km s -1 in (a), 
(b), and (c), respectively, and at 14. 4n 15 mjy beam -1 in (d). Crosses are at continuum peaks and filled ellipses show synthesized beams. 
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Fig. 2. — Molecular-line spectra at the two nuclei of Arp 220. Continuum has been subtracted. The right ordinate is the fractional 
absorption depth / a with respect to the continuum intensity of each nucleus, and can be converted to optical depth via r a = — log(— / a ). 
The la noise is 32, 12, and 28 mjy beam" 1 for HCO+ (4-3), (3-2), and CO(3-2), respectively. The CO(3-2) spectra overplotted in black 
are from 0'.'5-resolution data in Sakamoto ct al. (2008). The systemic velocity that we estimated for each nucleus is shown by a dashed 
line. 
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Fig. 3. — Position-velocity diagrams of Arp 220 nuclei. Position is measured from each nucleus along the position angle of 45° for the 
eastern nucleus and 270° for the western. Contours are at 2na (n = ±1,±2, ...) with negative ones dashed. The rms noise in brightness 
temperature is 3.1, 1.3, and 3.5 K for CO(3-2), HCO"'"(3-2), and (4—3), respectively. The horizontal dotted lines are at the systemic 
velocities of the nuclei, and the black rectangles show the spatial and velocity resolutions. 



TABLE 1 

Observed Parameters of Arp 220 



Parameter 


East 


West 


All 


/tot 


S , CO(3-2) [Jy km s" 1 ] 


250 


301 


(7.0 ±1.1) x 10 2 


0.23 ± 0.06 


S 'hCO+(4-3) [ J y km s_1 l 


19 


61 


(1.1 ±0.2) x 10 2 


1.0 ±0.3 


S HCO+(3-2) [ J y km s_1 ] 


35 


63 


(1.0 ±0.1) x 10 2 


0.64 ±0.13 


•50.87 mm [Jy] 


0.18 


0.49 


0.68 ± 0.10 


1.0 ±0.2 


Sl.lmm [Jy] 


0.07 


0.20 


0.27 ±0.03 


0.82 ±0.15 


CO (3-2) maxT b [K] 


58 (66) 


40 (47) 






HCO+(4-3) maxT b [K] 


16 (24) 


26 (34) 






HCO+(3-2) maxT b [K] 


12 (17) 


15 (21) 






0.87 mm maxT b [K] 


16 (23) 


42 (50) 






1.1 mm maxTb [K] 


7(12) 


20 (25) 






V svs (radio, LSR) [km s" 1 ] 


5415±15 


5355±15 







Note. — We conservatively estimate flux-calibration la errors to be 10% and 15% for 0.86 and 1.1 mm bands, respectively. Line 
fluxes are integrated over 4800—5900 km s _1 [5200—5900 km s -1 for HCO + (3— 2)] including both emission and absorption. Peak brightness 
temperatures are measured at or around each nucleus without correcting for beam dilution or missing flux. The first one is Rayleigh-Jeans 
T b and the second in parenthesis is Planck T b . 

1 Fraction of the flux (or flux density) recovered in our 0''3-resolution data with respect to the total flux of Arp 220. The total fluxes are 
106 ±23 Jy km s" 1 for HCO+(4-3) ijGreve et al.ll2009h . (3.0 ±0.6) X 10 3 Jy km s" 1 for CO(3-2) in the central 14" l|Wiedner et al.l[200^ . 
152 ± 26 Jy km s -1 for HCO + (3—2) (Th is work. SMA sub-compact array observations with the minimum projected baseline of 9 m), 
0.7 ±0.1 Jy for 0.87 mm continuum | |S08I and reference therein), and 0.33 ± 0.05 Jy for 1.1 mm continuum (This work, SMA sub-compact 
configuration). 



